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Abstract 

Non-equilibrium and equilibrium fluid systems differ due to the existence of long-range corre¬ 
lations in non-equilibrium that are not present in equilibrium, except at critical points. Here we 
examine fluctuations of the temperature, of the pressure tensor, and of the heat current in a fluid 
maintained in a non-equilibrium stationary state (NESS) with a fixed temperature gradient, a sys¬ 
tem where the non-equilibrium correlations are especially long ranged. For this particular NESS 
our results show that (1) The mean-squared fluctuations in non-equilibrium differ markedly in their 
system size scaling compared to their equilibrium counterparts and (2) There are large, nonlocal, 
correlations of the normal stress in this NESS. These terms provide important corrections to the 
fluctuating normal stress in linearized Landau-Lifshitz fluctuating hydrodynamics. 


In recent years fluctuations in fluids 
maintained in non-equilibrium steady states 
(NESS), and in non-equilibrium fluids in gen¬ 
eral, have attracted a large amount of atten¬ 
tion. Of particular interest are the studies of 


Evans, Cohen and Morris (ECM) 


and Searles Q, and of Gallavotti and Cohen 
3|. Later related work was done by Jarzynski 


Q. 


Evans 


^ and by Crooks |5]. Particularly in the ear¬ 
lier work, a focus was on non-equilibrium cur¬ 
rents and entropy production. For example, 
if Pxyr is the time average of the microscopic 
stress tensor, P^y, over a time interval r then 
in a NESS with a steady shear rate, 7 , ECM 
studied the probability distribution V{PxyT) 


in a particle system. Using an analysis 
based upon their computer simulations, they 
showed that, 

= drP^Mr)]] 

( 1 ) 

On the RHS the phase space time- 
dependence r{T) is restricted such that only 
the dynamical states i with a given value of 
Pxyr = ^ dTPxylPiir)] are taken into ac¬ 
count. They also numerically determined 
PiPxyr) and found it had a sharp peak 
around it’s average value and that although 
the bulk of the distribution was consis¬ 
tent with positive entropy production, there 
was a tail part consistent with negative en- 
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tropy production. The system size or N- 
dependence of the width of the distribution, 
V{PxyT) is an interesting open problem. 

Prior to this work, it was established in the 
1960's through the 1980's that correlations 
are very different in non-equilibrium fluid sys¬ 
tems than they are equilibrium systems 6|-l8|. 
For example, in the early 1980's it was pre¬ 


01 


dieted l9|l that the temperature and density 
correlations in a fluid in a NESS with a tem¬ 
perature gradient were extraordinarily long- 
ranged (in a sense growing with system size). 
This prediction was subsequently confirmed 
|ni . Ill with great precision in small angle 
light scattering experiments. For very small 
angle scattering, the scattering was found to 
be larger than the equilibrium scattering by 
a factor of 10^. All of this implies that the 
statistics of fluctuations in non-equilibrium 
fluids will in general be very different than 
those for fluctuations in the same fluid in an 
equilibrium state. For reviews see 


11 


13| 


Here we will expand on this point by ex¬ 
amining the system size dependence of not 
only the temperature fluctuations and distri¬ 
bution, but also the fluctuations and distribu¬ 
tions of the pressure tensor, and heat current 
for a fluid in a NESS with a temperature gra¬ 
dient. This particular NESS (see also point 6 
in the discussion) is unique because the cor¬ 
relations are so strong that they extend over 
the entire system size, see Eq. ([3]) . This is very 


different then the correlations in other NESS, 
which typically decay as a power law. We find 
that that temperature and pressure or nor¬ 
mal stress tensor fluctuations in a fluid with a 
temperature gradient are so large that many 
of the assumptions that are commonly made 
about non-equilibrium fluids when one uses 
equilibrium-like methods, such as maximum 
entropy formalisms, power series expansions, 
etc. to describe their properties, are called 
into question and may not be justified. We 
show that relative root mean square fluctua¬ 
tions of some thermodynamic and transport 
quantities do not in general scale with par¬ 
ticle number, as 1/N^P as do their equi¬ 
librium counterparts. We also describe the 
distributions of these fluctuations, and point 
out their anomalous properties as well. 

The new results in this paper are: (1) We 
show that the previously computed temper¬ 
ature fluctuations for a fluid in a tempera¬ 
ture gradient scale as and not as 

(2) We compute for the first time the 
pressure or stress fluctuations for a fluid in a 
temperature gradient. These fluctuations are 
shown to scale as (3) We compute 

for the first time the heat current fluctua¬ 
tions for a fluid in a temperature gradient. 
These fluctuations do not exhibit anamolous 
scaling, (4) We discuss the implications of 
these fluctuation results for the distributions 
of temperature, pressure, and heat current. 
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(5) We argue that the linearized about the 
NESS Landau and Lifshitz fluctuating hy- 
drodyanmics are not consistent with the long 
range nature of the normal stress fluctua¬ 
tions. 

We start with the well-known expression 
for the small wavenumber behavior of the 


temperature fluctuations 





ksT 
-|- Dr) 


Here p, v and Dt are the mass density, the 
kinematic viscosity, and thermal diffusivity of 
the fluid. This result is valid as long as k > 
l/L, with L the system size. Here we have 
assumed that the temperature gradient is in, 
say, the z-direction so that fc|| = -^/kl + ky is 
the magnitude of the wavenumber parallel to 
the confining walls. We note that this cor¬ 
relation function is long ranged as indicated 
by its k~‘^ behavior at small wave numbers, 
while as the equilibrium temperature fluctu¬ 
ations are localized in space with no singular 
behavior of the corresponding Fourier trans¬ 
forms at small wave numbers. In order to de¬ 
termine the spatial correlation of the temper¬ 
ature fluctuations one must invert the Fourier 
transform to obtain the correlation function 
Ct,ness{'<^) |lli- Using Eq. ([2]) we And that 
the spatial average of the temperature cor¬ 
relations has a strikingly different size de¬ 
pendence than the corresponding equilibrium 


temperature correlation. That is 

^T,NESS oc ——(3) 
pDriiy + Dt) 

where 


A 


T,NESS 


1 

1 ? 


drCT^NESsi'^), ( 4 ) 


where L is the characteristic system size. In 
general we will suppress numerical factors in 
spatially averaged quantities. The N depen¬ 
dence of this spatial average is easily found 
to be of order iV^^^'^^/'^while the equilibrium 
quantity would be 


At, eg — 


kuT^ 


( 5 ) 


L'^pCv ’ 

where is the specific heat per mass at 
constant volume. The A-dependence quoted 
above for Eq.([3]) is for fixed VT. Physically, 
and experimentally, it may be more reason¬ 
able to use VT = ^ with AT the tem¬ 
perature distance between the fluid confin¬ 
ing plates in the direction of the temperature 
gradient and fix AT. In this case the scaling 
of Eq.® is At,NESS oc oc i/N^d- 2 )/d_ 

In either case At,ness is large compared to 
At, eg in the large system size limit. Below 
we generally fix AT, to examine the system 
size dependence. The theoretical results ap¬ 
pear even more anomalous if the temperature 
gradient is fixed. 

For three dimensional systems a natural 
length, /, that occurs is 

ksT 


I = 


pDt{e Dt) 


( 6 ) 
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For water at STP, I ^ 3 x 10“®cm. The ra¬ 
tio of the non-equilibrium to equilibrium tem¬ 
perature fluctuations is 


^T,NESS 

^T,eq 



with a a molecular diameter. This ratio is 
generally large if L 3> a. Taking the reduced 
density to be unity, for water this requires, 
roughly. 


L ^ f 

a AT \ pcyi ) 

or, L/cr > 10 for typical experiments 
where AT/T < 1/5. This sets the scale of 
the system size where non-equilibrium fluctu¬ 
ations become dominant. If we assume that 
the temperature fluctuations in the NESS 
have a Gaussian distribution, we hud, again 
for three dimensions, that 



Next we consider pressure or normal stress 
fluctuations in this NESS. It is important 
to note that the average pressure or normal 
stress at some point in the fluid has a non¬ 
equilibrium term that is proportional to the 
square of the temperature gradient. This fol¬ 
lows from direct calculations, but it can be 
understood by symmetry arguments as well, 
since the pressure, for example, is a scalar 
quantity. The result for the most important 
non-equilibrium part of the (spatial and ther¬ 
mal) average pressure in three-dimensions is 
singular in the limit of large systems and is 
found to be fl Q 


P^Ess^AiL {VTf 



( 11 ) 


with 


VNESsi^T] ~ exp[-— (^) 

The corresponding probability distribution 
for equilibrium temperature fluctuations is 

~ ( 10 ) 

From a comparison of these two distributions 
one can see that the non-equilibrium distri¬ 
bution is dominant whenever the inequality 
given by Eq.(IH]) is satisfied. From Eq.(l9]) we 
see that in non-equilibrium the temperature 
fluctuations scale as 6T ~ ~ 

in three-dimensional systems (compared to 
l/jyi/s equilibrium systems). 


A = 


_ pcp(-y-l) 
2T 


1 / Ap\ I 1 /da) 

I I V ST I P 


dT J 


( 12 ) 


Here Cp, 7 , a are, respectively the specihc heat 
at constant pressure, the ratio of specific 
heats, and the coefficient of thermal expan¬ 
sion. More generally, the long distance part 
of the fluctuating pressure is |lJ, ll5| 


P(x) = A[6T{^f 


(13) 


Averaging this over the NESS and using 
Eq.(l 2 ]) gives Eq.lfTTl). 

The pressure fluctuations are 
<5P(x) = A[(5T(x))2 - {{5T{x)f)EEss] 
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and the non-equilibrium correlation function 
is Cpp(x,y) =< (5P(x)5P(y) >ness- Using 
a Gaussian approximation, we find this 
correlation function to be, 

Gpp(x,y) = 2/l2[< 5T(x)5T(y) 

>nessT 

(14) 

The growth of the temperature correlations 
with distance implies that the spatially aver¬ 
aged Cpp behaves as, 

^P,NESS OC [- — -] OC — (15) 


fluctuations scale as SP 1/L l/iVU3 
(compared to in equilibrium systems). 

As a final example we consider the NESS 
fluctuations in the heat current in some di¬ 
rection, 6jH,a, where a = {x,y,z). There 
are both linear and nonlinear (analogous to 
Eq. (lT3l) ) contributions to the heat flux. The 
linear contribution can be obtained by identi¬ 
fying the microscopic fluctuating heat current 
as 5jp'(x) = —AV5T(x) with 6T the fluctu¬ 
ating temperature gives. 


while in equilibrium, 

nkpT (dP\ 1 


Again, the non-equilibrium fluctuations are 
large compared to the equilibrium ones at 
large length scales. Comparing Eqs. fllSp and 
(USD and using Eq.(l8]), this roughly occurs 
when Lja > e^, or, in typical experiments, 
Lja > lOh 

If we assume Gaussian behavior then the 
normal stress fluctuation distributions in the 
NESS and in equilibrium are 


< ST(x)ST(y) >HESS 

(19) 

The long ranged part of the temperature fluc¬ 
tuation in Eq.([T9]) is given by Eq.(l2]), with the 
final result for in Fourier space given by. 

This result can also be directly obtained with 


/ §p 

VxEssm ~ exp[-- 


(17) 


and 

Veq{SP) ~ exp 


(A] ^spf 


2nkBT \dP 


(18) 


That is, in three-dimensional non¬ 
equilibrium systems the normal stress 


kinetic theory methods 1^. Here the k de¬ 
notes unit vectors. The dependence in 
Eq. (l20|) implies that for systems in d dimen¬ 
sions, —)■ oo) ~ that is, power 

law correlations. The angular factors don’t 
really change this in general. For example, it 
is easy to see that if in three-dimensions one 
look in the middle of the fluid at = Uz and 
considers correlations along that plane (per¬ 
pendicular to the direction of the tempera¬ 
ture gradient) as a function of rp one finds 
that the z—component of the current fluctu¬ 
ations is ~ l/r^. Physically this means that 
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if there is a current fluctuation in the wrong 
direction then the fluctuations in that entire 
plane will likely also be in the wrong direc¬ 
tion. To obtain a better measure of these 
correlations we consider the spatial average 
of Ci^(x), 


A 


kfiT c„X 


{ATf 


^ CpA 2-d(yrr..2 ^ 

H.NESS ^ ^ (VJ) OC 

( 21 ) 

That is, even though the current correlations 
are of long-range, they are weighted by a fac¬ 
tor of (VT)^ cx 1/L^ , for hxed AT, so that 
qualitatively scales just like Ajj^eq l/T*^ 
and is, in fact, small compared to Ah, eg be¬ 
cause it is of relative order ■ 

In analogy with Eq. ifT^ . the non-linear 
portion of the fluctuating heat current is 
jni^) oc 5T(x)u(x), with u the fluctuating 
fluid velocity. The fluctuation correlation of 
this quantity gives C'^'^(x). It can be read¬ 
ily computed (ol, luj and in space it decays 
as ~ ~ {ATf /[L‘^r‘^^-% To 

obtain a better measure of these correlations 


we consider the spatial average of C 


NL( 
H \ 


nr 

^H,NESS ^ 


(AT)2 

l2d-2 


( 22 ) 


These correlations are even weaker then those 
from the linear portion of the fluctuating heat 
current and can therefore be neglected. 

We conclude with a number of remarks: 


1. The result for Aptness, Eq. ifT^ . has 
an important implication. In Landau 


and Lifshitz (L&L) fluctuating hydro¬ 
dynamics an input is that the fluc¬ 
tuating stress tensor is local in space 
and time. Averaging the normal com¬ 
ponent of the L&L fluctuating stress 
tensor correlation function in three- 
dimensions over space and and a mi¬ 
croscopic time interval leads to a quan¬ 
tity A that scales like A cx 1/L^ which 
should be compared with Aptness oc 
1/T^. That is, our results indicate that 
there are important corrections to the 
usual linearized L&L equations due to 
the very long-range correlations that 
exist, in at least some NESS. 

In a non-linear fluctuating hydro- 
dynamic description these anomalous 
terms will be generated by fluctuation 
or renormalization effects. 


2 . Compared to the pressure or stress fluc¬ 
tuations both the linear and non-linear 
heat current fluctuations are of rela¬ 
tively short range. This result is ac¬ 
tually important: It is consistent with 
the local assumption for the fluctu¬ 
ating heat current in the L&L ap¬ 
proach. This in part explains why con¬ 
ventional linearized fluctuating hydro¬ 
dynamics can be used to successfully 
compute the very long range NESS 
temperature fluctuations 111. llTf. 
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3. As another measure of just how anoma¬ 
lous the normal stress fluctuations are, 
note that if a fluctuation as large as the 
average value of Pness is considered, 
Eq. flTT]) . in Pness[^P] fhen the prob¬ 
ability of that fluctuation is of 0{L^). 
That is, the normal stress is not a self¬ 
averaging quantity. 


4. The correlations in a fluid with a con¬ 
stant velocity gradient (that is, the 
case studied by ECM) are very differ¬ 
ent (see, for example, ISj) then either a 
fluid with a temperature gradient, or in 
a mixture with a concentration gradi¬ 
ent. Although there are power law cor¬ 
relations in the velocity gradient case, 
they are sufficiently weak that the nor¬ 
malized NESS stress-tensor correlation 
function does scale like 1/NP^. This 
is true whether the velocity gradient, 
Vm, with u the hydrodynamic velocity, 
is fixed or if Vu = Au/L is used, and 
Au is fixed. 


5. In our presentation we have ignored a 
discussion of the boundaries, that must 
be present in this NESS. For the case of 
perfectly conducting walls the temper¬ 
ature fluctuation problem in a bounded 
geometry has been discussed and solved 
elsewhere jll|. 


given here for a single component fluid 
with a temperature gradient are ob¬ 
tained in a fluid mixture with either 
a concentration gradient or a temper¬ 
ature gradient (see jn| and references 
therein). In this case, the temperature 
fluctuations, the concentration fluctua¬ 
tions, and the normal stress, or pres¬ 
sure, fluctuations are all found to obey 
anomalous statistics. 


7. The profound difference between cor¬ 
relations in non-equilibrium and equi¬ 
librium systems means that many of 
the theoretical techniques developed 
to describe equilibrium systems cannot 
be applied to non-equilibrium systems. 
Non-equilibrium quantities do not have 
virial expansions g. A local expan¬ 
sion of the fluxes or currents in terms 
:he gradients is also not 


of powers c 

)f t 

le g 

possible js. 

19, 

201 


20|. Other techniques 


such as maximizing an entropy (for ex¬ 
ample, the so-called max cal method 


2 lj) to obtain a non-equilibrium distri¬ 


bution function may not work, at least 
in their most naive form. 
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6 . Structurally identical results as those ber DMR-1401449. 
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